Number Fluctuations and Energy Dissipation in Sodium Spinor Condensates 
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We characterize fluctuations in atom number and spin populations in _F = 1 sodium spinor 
condensates. We find that the fluctuations enable a quantitative measure of energy dissipation in 
the condensate. The time evolution of the population fluctuations shows a maximum. We interpret 
this as evidence of a dissipation-driven separatrix crossing in phase space. For a given initial state, 
the critical time to the separatrix crossing is found to depend exponentially on the magnetic field 
and linearly on condensate density. This crossing is confirmed by tracking the energy of the spinor 
condensate as well as by Faraday rotation spectroscopy. We also introduce a phenomenological 
model that describes the observed dissipation with a single coefficient. 

PACS numbers: 03.75.Mn, 67.85.-d, 03.75.Kk 
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The transition from a thermal atomic gas to a Bose- 
Einstein condensate (BEG) is marked by the appearance 
of a scalar order parameter. Spinor BEGs have an ad- 
ditional spin degree of freedom which results in a vector 
order parameter. The increase in complexity leads to the 
formation of spin domains [il, B 3, the appearance of 
novel phases [J] and the possibility of high spatial reso- 
lution magnetometry Spin-1 spinor BEGs have been 
studied with ^'^Na atoms that show antiferromagnetic in- 
teractions 0, @, 0, Q and with *''Rb atoms that show 
ferromagnetic interactions 0, [l^, 11, 12, [l^l- A re- 
markable result is the observation of spin population os- 
cillations that appear when the system is taken out of 
equilibrium in the presence of a magnetic field [1,01. An 
interplay between the quadratic Zeeman energy and a 
spin-dependent interaction energy determines the oscil- 
lation frequency. The oscillations are nearly harmonic 
except near a separatrix in phase space where the period 
diverges 0,i,[iai- The system can be forced onto either 
the low or high energy side of the separatrix using the 
magnetic field strength 0, [l3| , the BEG density Q , or 
the balance among spin states Q . 

Quantum optical effects in spinor BEGs arc now being 
actively studied. Recent experiments observed the spin- 
mixing analogue of parametric amplification IJ] . Num- 



ber fluctuations and spin oscillations were investigated 
in a ferromagnetic Rb BEG in Refs. and 18. The ob- 
served spin oscillations damped and the system reached 
a steady state, while the fluctuations saturated. Non- 
dissipative theories of quantum effects in number fluctu- 
ations [3, [13, [3] show that such damped spin os- 



cillations can be produced by dephasing from quantum 
fluctuations or from number and phase fluctuations in 
the initial state. 

We report observations of the dynamics of atom num- 
ber fluctuations in an antiferromagnetic ^"^Na spinor BEG 
which show strong evidence of energy dissipation. Atom 
number fluctuations of the spin projections are extracted 



from a scries of measurements on condensates which spin- 
mix while slowly evolving to the ground state. Because 
of the dissipation, a spinor BEG with given initial con- 
ditions and sufflciently high energy will cross the sep- 
aratrix at a critical time, tc- The time evolution of 
the population fluctuations unambiguously identifles tc- 
We have developed a dissipation model that includes 
a single phenomenological coefficient, is classical, and 
docs not require the intrinsic quantum fluctuations that 
have been used to describe ferromagnetic Rb spinor sys- 
tems. Other dissipation mechanisms have been suggested 
m but do not explain our data. Mean-fleld 

simulations at zero mag netic field including finite temper- 
ature effects in Ref. |2l| indicate that thermal excitations 
play a prominent role in the dissipation. 

We show that a non-dissipativc quantum model does 
not account for our data. This model predicts a damped 
spin oscillation due to quantum dephasing similar to the 
ferromagnetic case. Its predicted steady states, however, 
are very different from the experimental observations in 
this paper and Ref. @. It is interesting to note that a 
spinor BEG represents a nominally isolated quantum sys- 
tem that shows dissipation and does not display quantum 
rephasing over our observation times. 

The setup is similar to that of our previous work 0]. 
We create a F = 1 spinor BEG of TV = 1.50(3) x 10^ '^'^^a 
atoms through 6 s forced evaporation in a crossed optical 
dipole trap using a multi-mode fiber laser at 1070 nm (all 
quoted uncertainties are one standard deviation, com- 
bined statistical and systematic) . Wc apply a weak mag- 
netic field gradient during the evaporation to fully po- 
larize the atoms to the \F = l^mp = +1) state. The 
final trap oscillation frequencies are wqCv^, 1, 1) in the 
three spatial directions. We have performed experiments 
with measured values of (x'o/(27r) = 154(5) Hz, 220(7) Hz, 
and 305(9) Hz. This corresponds to mean Thomas-Fermi 
radii of 6 /xm to 8 /xm. The density of the BEG is changed 
only by changing the trap frequency and not by reducing 
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the atom number, as this causes a large loss in signal. 

To prepare the initial state, we turn off the magnetic 
field gradient and ramp to a magnetic field, B, less than 
61 /iT. We apply an rf pulse resonant with the linear Zee- 
man splitting (frequencies of hundreds of kHz) to rotate 
the atomic spin. All of our experiments start with the 
same initial atomic state which has p+i = p-i = 1/4 
and po = 1/2, where the pmp are fractional populations 
for the three Zeeman sublevels {mp = 0, ±1). This state 
has zero magnetization, where the magnetization is de- 
fined as m = p+i — p_i. We use two methods to detect 
spin mixing dynamics: Faraday rotation spectroscopy 
and Stern-Gcrlach separation combined with absorption 
imaging (SG-AI). Faraday rotation spectroscopy can be 
used for continuous observation of spin oscillations of a 
single BEG over short time scales, while SG-AI can di- 
rectly measure spin populations, albeit destructively Q. 

The Mandel Q parameter is a common way to quantify 
fluctuations in quantum systems, with Q > (< 0) rep- 
resenting super- (sub-) Poissonian distributions [24|. We 
use a modified Q parameter to characterize population 
fluctuations during spin mixing. This Q parameter of pQ 
is dcflncd as llSll 
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where (N) is the mean value of the atom number in the 
BEG. (Apg) and (po) correspond to the variance and 
the mean value of po, respectively. At each delay time 
after initialization, we extract the variance from 25-30 
repeated SG-AI measurements. We use po rather than 
the population, because measurements of po and m arc 
less sensitive to the 2% fluctuations in the initial N. 

Theoretically, our initial state, prepared from a single 
component BEG, should be a coherent state with a Pois- 
sonian atom number distribution (Q — 0). The observed 
Q at < = is equal to 4, while the minimum observable Q 
at other times depends on the populations, due to tech- 
nical noise in atom counting. In particular, when one of 
the spin populations is close to zero the minimum value 
of, and the error in, Q are larger. 

Figure [1] shows the time evolution of Q for m = at 
four magnetic fleld strengths. Q has a value equal to the 
experimental limit at t = 0, increases to a peak at tc, 
then decreases back to the experimental limit within 10 
s. We find that at t = 10 s all the remaining atoms are 
in the mp = state (~20% are lost). A Gaussian fit to 
Q{t) is applied to extract tc- We observe that the value 
of tc decreases with increasing field. For high fields (e.g., 
B = 60.7 fiT in Fig. [T|) where the system is on the low 
energy side of the separatrix at t = 0, Q{t) stays at the 
experimental limit during the whole evolution. 

The single mode approximation (SMA) [l3l appears to 
be a suitable model to explain our data. In this approx- 
imation all the spin components share the same spatial 
wavefunction <f>(r) and the total wavefunction is '^{r, t) — 
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FIG. 1: (Color online) Time evolution of Q with (n) = 
1.37(6) X 10" cm"^ {c/h = 33(1) Hz) at four magnetic fields 
as indicated in the figure. Inset: Time evolution of i?iimit /Ssep 
for the same data (see text). Good agreement between the 
predictions of Eqs. (3) and the data is found using a /3 that 
falls on the line shown in the inset of Fig. [3] An example 
curve is shown for B = 21.2 pT (thick black line). 



$(r)(y^I7(i)e»«-(*), yMi)e''"^*^ ^/^^)e^'^-(')), 
where Omp represents the phase of each spin component. 
Taking into account the conservation of m and 
the description simplifies into a model with only two 
dynamical variables po and 9, where 9 = 9^i+ 9^i — 29q. 
The classical spinor energy is [l^l 



E = E^,(l- 



{l-po)+cpo ({1 - Po) + V(l - Po)'^ - m'^ cos 6*^ , 

(2) 

where Eq^ oc B^ is the quadratic Zeeman shift [Eq^/h = 
(0.0277Hz/(^T)2)B2), c = C2{n) is the spin-dependent 
interaction energy with the mean BEG density (n), and 
C2/h = 2.4 X 10^13 Hz cm^ for 23Na {h is the Planck 
constant) Q. In the Thomas-Fermi approximation, (n) 
and c are proportional to N'^^^lUq^^. The separatrix is 
the contour in (po, 9) phase space with energy i?scp, on 
which there is a saddle point where p'o = ^ = 0. For an 
antiferromagnctic spinor BEG with m = 0, Egep = Eq^. 
The mean- field ground state is po = 1 for m = 0. 

The energy of the system cannot be directly inferred 
from SG-AI measurements. We can, however, use the 
values of po and m to calculate an upper bound (i?iimit) 
to the classical spinor energy. In fact, -Bumit is equal to 
(1 - Po'''')[Eq^ + 2cpg^^'^) for TO = 0, where pg^'^'' is the 
maximum value of po among the 25-30 repeated measure- 
ments. The inset of Fig. [T] shows that iJumit decreases 
over time and crosses the separatrix around tc at low 
fields. This indicates that E is not conserved. 

One source of dissipation is loss of BEG atoms and the 
corresponding decrease in density, resulting in a decrease 
in c. The interaction energy evolves as c{t) oc e"^'''*/^ as 
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FIG. 2: (Color online) Time evolution of Q with in) = 
2.02(8) X 10" cm^^ {c/h = 48(2) Hz) aX B = 21.2(1) pT. 
Red dots represent experimental data. The dashed (black) 
line is the result of the non-dissipative quantum simulation. 
The solid (blue) line represents a classical Monte Carlo simu- 
lation using Eqs. (3) with /3 7^ 0. Inset: A classical simulated 
path including energy dissipation through phase space. The 
black dot represents the initial state. The (blue) dashed line 
is the contour with energy iJsep. 



follows from the Thomas-Fermi approximation, with 7 < 
0.02 s~^ estimated from the observed atom loss. Because 
all observed tc are much smaller than I/7, the evolution 
of c does not explain our data and, hence, is ignored. 

In Fig. [5] the evolution of Q is compared to several 
theoretical models. The dashed line is a result of a non- 
dissipative quantum simulation based on a quantized ver- 
sion of Eq. [2] following the prescription in Ref . |l8|. An 
initial Gaussian wavepacket with a standard deviation 
in po of 0.8% mimics the fluctuations in initial popula- 
tion. A classical Monte Carlo simulation based on Eq. [2] 
provides a similar result (not shown). We average over 
30 trajectories using a Gaussian probability distribution 
for the same variance of po- These two non-dissipative 
simulations approach an identical steady-state value of 
Q, however, they do not explain the observed Q. More- 
over, the steady-state spin populations of these models 
are different from experimental observations. 

We modify the equations of motion for po and 6 
by adding a dissipation term inspired by the description 
of ohmic loss in Joscphson junctions [23|. This leads to 



po = -{2/h)dE/de, 

9 = +{2/h)dE/dpo + f3po, 



(3) 



which corresponds to an evolution of the energy given by 
dE/dt = -h(3{po)^/2. The solid line in Fig. [g] represents 




FIG. 3: (Color online) Measured tc as a function of 
B at mean densities, 2.02(8) xlO" cm"^ (blue circles), 
1.37(6) X 10" cm"^ (red squares) and 0.89(4) xlO" cm"^ 
(black triangles). The solid lines are exponential fits. Inset: 
/3 as a function of -Bqz/c for the same data. The solid line is a 
weighted linear fit to the data. At the separatrix Eqz/c = 1. 



a classical Monte Carlo simulation with this dissipation 
mechanism. We sample two Gaussian distributions with 
standard deviations of 0.8% and 1% for the initial values 
of po and c, respectively. The deviation in c is due to 
the 2% uncertainty in N, while the drift in luq over these 
measurements is negligible. 

By assuming E(tc) = Esep the coefficient /3 is obtained, 
where tc is the observed critical time. For an extracted /?, 
we find agreement between the observed -Enmit and the 
energy derived from the dissipative model, as shown in 
the inset of Fig. [TJ Additionally, including the initial fluc- 
tuations enables our model to reproduce the behavior of 
Q{t) in Fig.O The predicted Q, however, is almost twice 
as large as the experimental value at tc- Several other 
phenomenological dissipation terms were tested, but the 
term (3po is the only one linear in po, 9, p'o, or 9 that we 
found to drive the system to the correct ground state 0] . 

We can understand the evolution of Q from consider- 
ing the trajectory in the inset of Fig. [2l It shows a sim- 
ulated path through phase space based on Eqs. [H The 
dissipation leads to a gradual decrease of energy, which 
results in a larger oscillation amplitude in po as the tra- 
jectory approaches the separatrix. Paths starting from 
slightly different initial conditions separate after several 
oscillations. Larger oscillation amplitudes in po in these 
separated paths lead to a larger variance and thus an in- 
creasing Q. As the energy decreases past the separatrix, 
the oscillation amplitude decreases until it goes to zero 
at the ground state. This reduces Q back to a minimum. 
Thus Q reaches a maximum as the ensemble crosses the 
separatrix at time tc- The energy dissipation is determin- 
istic and does not explain the observed Q by itself. The 
evolution of Q is a consequence of dcphasing due to the 
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FIG. 4: Faraday signal for a single evolving BEG with (n) = 
1.37(6) X 10^'' cm"^ (c/h = 33(1) Hz). For this specific case, 
the initial po is 0.45. The transition from the oscillating phase 
solution to the running phase solution happens between 35 ms 
and 45 ms, as shown by the signal periodically approaching 
zero at the minima afterwards. This is in good agreement 
with the predicted tc from Fig. O 



spread of initial conditions and the change in oscillation 
amplitude during the evolution of the system. 

Figure [3] shows the measured tc as a function of B for 
three densities (and thus for three values of c). The solid 
lines show exponential fits to the data. Uncertainties are 
larger for data at the smallest density because we find 
that Q{t) shows a broader and asymmetric peak and thus 
it is hard to extract tc from a single Gaussian fit to Q{t). 

The inset of Fig. [3] shows /3 determined from Eqs.[3]as 
a function of E^z/c. The coefficient (3 fits to a linear func- 
tion with a (dimensionless) slope of 0.22(1). The relation 
holds for different B and (n), hence, a single function 
seems to describe all the data in terms of a single param- 
eter Eqz/c. It is important to note that this same ratio 
completely determines the shape of the energy surface in 
phase space for a given to. For very small fields, the fit 
extrapolates to an unphysical negative (3. This indicates 
that either the dissipation term or the functional form 
of (3 are perhaps incomplete or not appropriate. For ex- 
ample, our analysis does not include the coupling of the 
single BEG mode to other degrees of freedom, such as 
elementary excitations (whose energies are on the order 
of hojo), which are indicated to be important by Ref. |21. 

We have also studied the phase-space dynamics using 
Faraday rotation spectroscopy. This method can reveal 
evidence of the separatrix crossing from a single BEG 
realization. As outlined in Ref. d, a BEG on the high en- 
ergy side of the separatrix produces an oscillating Fara- 
day signal with non-zero minima, while on the low en- 
ergy side the Faraday signal reaches zero. Figure |4] shows 
an example of the Faraday signal from a single BEG at 
B = 36.3(1) yuT and m = 0. The trace shows a tran- 
sition from non-zero minima to minima near zero, thus 
providing a clear signature of crossing the separatrix in 



phase space between 35 ms and 45 ms. While the de- 
tails of repeated traces vary, the system always crosses 
the separatrix during this narrow time interval. 

Figure 5] confirms our assumption that tc corresponds 
to a crossing of the separatrix, and tc from the Faraday 
signal agrees with the extrapolated value of from Fig.[3l 
Light-induced atom loss limits our Faraday detection to 
100 ms, making it hard to observe the separatrix crossing 
using Faraday detection for smaller values of B. 

In Fig. |4] the rapid reduction of the oscillation ampli- 
tude with time is due to atom losses generated by the 
Faraday beam and not, as might be expected, to energy 
dissipation. The Faraday beam gradually destroys the 
spin dynamics by off-resonant light scattering and tensor 
light-shift dephasing, and thus the Faraday detection is 
only effective over short observation times. 

In conclusion, we have studied spin population fluctu- 
ations and energy dissipation in a spinor BEG. Popula- 
tion fluctuations peak at a critical time where the en- 
ergy of the system equals that of the separatrix in phase 
space and we have confirmed the separatrix crossing us- 
ing Faraday rotation spectroscopy. We present a dissi- 
pation model with a single phenomenological coefficient 
that describes our data. While the underlying physics 
requires further study, this work sheds new light on dis- 
sipation mechanisms in a nominally-isolated spinor BEG. 
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